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1 Introduction 

Typically, dark energy (DE) models are based on scalar fields minimally coupled to gravity, 
and do not implement the explicit coupling of the field to the background matter [1, 2]. How¬ 
ever there is no fundamental reason for this assumption in absence of underlying symmetry 
which would suppress the coupling. Given that we do not know the true nature of DE and 
dark matter (DM), one cannot exclude that there exists a coupling between them. Whereas 
new forces between DE and baryon matter particles are heavily constrained by observations 
(e.g. in the Solar system and gravitational experiments on Earth), this is not the case for DM 
particles. In other words, it is possible that the dark components interact with each other, 
while not being coupled to standard model particles. The study of the interaction of DE 
and DM is an important and promising research direction [3-5]. Moreover, disregarding the 
potential existence of an interaction between dark components may result in misinterpreta¬ 
tions of observational data [6]. Since there is no fundamental theoretical approach that may 
specify the functional form of the coupling between DE and DM, presently coupling models 
are necessarily phenomenological. Of course, one can always provide arguments in favor of 
a certain type of interaction. However, until the creation of a microscopic theory of dark 
components, the effectiveness of any phenomenological model will be defined only by how 
well it corresponds to observations. 

Interaction between the dark components is phenomenologically described by the fol¬ 
lowing modification to the conservation equations 

Pdm T ‘^Hpdm Qi Pde 3i4(l -|- w'jpde (l'^) 

where H is Hubble parameter, pdm and p^e are densities of DM and DE components respec¬ 
tively. The interaction function Q depends on the scale factor and can take any possible 
form Q = Q {H, Pdm, Pde^t) ■ However, physically, it makes more sense that the coupling be 
time-independent with preference given to the factorized-i4 dependence Q = 2>Hq{pdm, Pde)- 
We consider the simplest linear models with the interaction functions 

Ql Pdm, Qll Pde, QlII (^Pdm T Pde), (1'^) 

which allow exact solutions discussed below, where we find model-independent^ constraints 
on the coupling parameter 5 between DE and DM in the models (1.2) formulated exclusively 
in terms of the cosmographic parameters (CP) values, which can be directly extracted from 
the cosmological observations. 

'^in the sense of independence w.r.t. choice of particular gravity model 
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2 Cosmography background 


Cosmography is an approach entirely based on the cosmological principle, stating that the 
Universe is homogeneous and isotropic on scales larger than a hundred megaparsecs. It 
allows to choose among whole possible variety of models describing the Universe a narrow 
set of homogeneous and isotropic models. The cosmological principle enables us to build 
the metrics and make first steps towards interpretation of the cosmological observations. 
Cosmography is just the kinematics of cosmological expansion. In order to build the key 
characteristic — time dependence of the scale factor a{t) — one needs to take equations of 
motion (the Einstein’s equations) and make an assumption about material content of the 
Universe, which allow to construct the energy-momentum tensor. Cosmography is efficient 
because it allows to test any cosmological model which does not contradict the cosmological 
principle. Modifications of General Relativity or introduction of new components (such as 
DM and DE) evidently change the dependence a{t) but it does not affect relations between 
the kinematic characteristics. 

In order to make more detailed description of kinematics of cosmological expansion it is 
useful to consider the extended set of the parameters which includes the Hubble parameter 
H{t) = a ja and higher order time derivatives of the scale factor [7-9] q{t) = —C 2 , j{t) = C 3 , 
s{t) = 6 * 4 , l(t) = Cs, where 

C = 

” a df^ a dt 

Dunajski and Gibbons [10] proposed an original way to test the cosmological models based 
on the General Relativity or its modifications. The suggested procedure is the following: 

• transform the Friedmann equation to ODE for the scale factor; in order to do that, use 
the conservation equations for each component included in the model to find explicit 
dependencies of the energy densities on the scale factor; 

• the obtained ODE includes a row of free parameters (such as initial values of the energy 
densities and the curvature); let their number equals N] differentiate this equation 
w.r.t. time N times and express higher time derivatives of the scale factor in terms of 
the cosmographic parameters (CP); 

• solve the obtained system of N equations w.r.t. the N free parameters and so express 
them through the CP; 

• substitute the obtained expressions for all N free parameters into the initial Friedmann 
equation to obtain the relation between the CP corresponding to the model under con¬ 
sideration; the precision which the obtained relation holds with will determine relevance 
of the considered model. 

Such approach can be easily generalized to the models which include such nontrivial effects 
as interaction between the components, volume and shear viscosity, decay of vacuum energy. 
Higher CP are presently known with insufficient precision, however perspective of their cor¬ 
rection in near future makes the proposed method an efficient tool for testing of cosmological 
models. 
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3 Cosmographic constraints on the coupling constant 

We consider a general case of Universe filled with two components labeled 1 and 2 with the 
EoS and the interaction function respectively 


Pi=wipi, P2 = W2P2, Q = 36Hpi. (3.1) 

In particular, the case Qj in (1.2) correspond to wi = 0, W 2 = w and the case Qu is obtained 
with wi = w, W 2 = 0 and 6 — —6. 

Transformation of the independent variable 


d _ 

— = aH— 
at da 

allows to exclude the time variable from (1.1) and rewrite it in the form 

ap'i + 3{l + wi - 6)pi = 0, ap2 + 3{1 + W2)p2 =-36pi, (3.2) 

where the prime denotes derivation w.r.t. the scale factor. The system (3.2) has exact solution 


Poi 


Pi(«)= 3(1+^!-^) ’ ^2(a) = 


P02 + dpol/{W2 - Wi + (5) 5poi/{W2 - Wi + (5) 


Q^o\^i.-rLui—uj - ■ q3{1-\-W2) q3{1-\-wi—5) 

where pQi = pi{a = 1), i = 1,2. Then the first Friedmann equation becomes 

■ 2 , ^ B 

aP 

where 87rG/3 = 1, a = 3(t(;i — 6) + 1, (3 = 3w2 + 1, 


^ = Poi 1 - 


1 


1 + 


B = P02 + 


dpol 


102—101 


W2 — Wi + 6 


Then we differentiate the equation (3.3) two times w.r.t. time to obtain 


—2— = 2 H^q = Aaa 


+ BPa 


-P-2. 


-P-2 


2 — = 2H^j = Aa(a + l)a"""^ + B/3(p + l)a 
a 

and solve the obtained system of two linear equations w.r.t. A and B 

t2(a (R I Ot 72 


_ (j - (/3 +l)g) «+2. „ ^ 2g [j -{a + l)q) 0+2 

a(a - /3) ’ /3(/3 - a) 


a{a — p) 

After substitution of (3.5) into (3.3) one obtains 


(3.3) 


(3.4) 


(3.5) 


^ = + ^ + (3,6) 

The third time derivative of the expression (3.3) leads to 

2 fa^^ I _ Aa{a + l){a + 2) i?/3(/3 + l)(/3 + 2) 

H\d ^2 ) {s + qj)- 
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(3.7) 


Using the expressions (3.5), the latter equation can be transformed to the form 

■s + + (« + /3 + 3)i - q{a + l)(/3 + 1) = 0. 

Substituting the explicit expressions for a and /3, (3.6) and (3.7) take on the form 

k _ 2[^q{w2 + - i5 + 1) - j] 

a^H'^ (3u; 2 + l)[3(rci — (i) + 1] 

s + qj + [3(u; 2 + rci - 5) + 5]j - q{2>W2 + 2)[3(r(;i - 5) + 2] = 0. 
From (3.9) one obtains 6 as explicit function of CP and EoS parameters 

^ 0,0 , s + qj+ ‘i{w2 + l)j 

3 - {3w 2 + 2)q 

Under the assumption of flatness (A: = 0), (3.8) gives an alternative relation 

^ j - q{3w2 + 2) 

3 ^ ^ 3[(3u; 2 + l)/2 - q] ' 

For the models Qj and Q// in (1.2) one obtains respectively 


A _ ^ I ^ 2 ~ 9(3^^ + 2) 

^ 3^ 3 3w + l-2q ' 

5/7- 3+^+3Y3^, 


51 — 2 + 


s + qj + 3{w + l)j 


Sjj — 2 3w -f- 


j - {3w + 2)q 
s + {q + 3)j 


j - 2g 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 


4 Test of the models with the observational data 

We plot 6iji and 5jjj (see figures 1 and 2) for —1.1 < w < —0.7 (it is the present ob¬ 
servational estimate on the EoS parameter for DE), taking the CP values Hq = 74.22^5 Qg, 
<70 = “0-61491 q 2238 ! Jo = 1-03015'qq^, '^o = 0.161}'q| from [11], where the authors use the 
data of the SNela Union 2 compilation by the supernovae cosmology project [12]. 

The model Qm also has exact but cumbersome solution: 


Pdm — 


Pde 


PdmO 

PdeO 


1 


^a+2 


tCk+2 


+ 


1 


+ 


1 


+ 


a 


0+2 


where 


{a, /?} = 1 -|- 3(5 


PdeO 215 ) PdmO / 

^ 1 


(4.1) 

2v/i + (i + 5)' 


a/5+2 J ’ 

PdmO PdeO / 

^ 1 


(4.2) 

2^1 + {! + §)" 


a/5+2 J ’ 

w , 1 w\'^ 





The corresponding Friedmann equation again has form (3.3), the relations (3.6) and (3.7) 
hold and give quadratic equations for 6, so one finally obtains (see figure 3) 



1 - 


9qw^ 


[s - 2q{2 + 3w) + j(5 + q + 3tc)] 


1 - 


9w'^ 


[1 -|- 2j + 3w — 6g(l + re)] 


(4.3) 

(4.4) 
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Figure 1. The EoS parameter dependence of the coupling constants S and S in the model Qj (3.12). 
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Figure 2. The EoS parameter dependence of the coupling constants S and S in the model Qn (3.13). 

Careful analysis of the plots reveals a remarkable feature: the functions 6{w) corresponding 
to different models {I, II, III) have exact common roots, namely: 

57 ( 1 (^ 01 ) = 5ii{wm) = 5i(rcoi) = 0, wqi = - ^ ~ -1.06; (4.5) 

3(d - 2q) 

— — — 1 -|- 2 ? — Qq 

Si{wo2) = Sn{wo2) = Si{wo2) = 0 , W02 = — 3(1 _ 2q) ^ (4.6) 

For the case of DE in form of cosmological constant {w = —1) the corresponding estimates 
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Figure 3. The EoS parameter dependence of the conpling constants 5 (4.3) and 5 (4.4) in the model 

Qiii- 


for the coupling constant 6 are the following: 

(5/(-l) Ri 0.86, (5//(-l) Ri-0.16, (5i(-l) «-0.06; (4.7) 

6 i{-l) « -0.026, Sn{-1) « -0.00897, 5i(-l) ~ -0.0066. (4.8) 

We disregard the solutions 62 and 82 from (4.3) and (4.4) respectively as they give unreason¬ 
able results for the coupling constant. 



5 Summary 

We considered three cosmological models with linear type of interaction between the dark 
components and obtained two types of restrictions on the coupling constant 6 in terms of the 
CP: the first type contains the CP q, j and s, the second is obtained in assumption that the 
Universe is flat {k = 0) and contains only q and j. We resolved the obtained constraints to 
explicitly obtain (5 as a function of CP and EoS parameter w and plotted the dependence in 
the interval —1.1 < w < —0.7 which corresponds to the present observations. We see that all 
the three considered models give similar results for the coupling constant 6 , with exceptions 
of the model Qj where 6 {w) diverges in the considered interval (for w ~ —0.74). Numerical 
estimates of the coupling constant for the case of cosmological constant {w = —1) for all the 
three considered models give small negative value of 6 , with exceptions of the model Qj where 
1 . 
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